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Abstract 

We study an extended Skyrme model which includes fourth and sixth-order terms. 
We explore some static properties like the A-nucleon mass splitting and investigate the 
Skyrmion breathing mode in the framework of the linear response theory. We find that 
the monopole response function has a pronounced peak located at ~ 400 MeV, which we 
identify to the Roper resonance iV(1440). As compared to the standard one, the extended 
Skyrme model provides a more accurate description of baryon properties. 
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INTRODUCTION 

At low energy the QCD running coupling constant a s becomes large which renders the 
standard perturbation theory inapplicable. In order to describe hadronic physics at low energy, 
effective theories including the main features of QCD (e.g. chiral symmetry) have therefore 
been proposed. These theories carry the label "effective" because their degrees of freedom are 
the hadronic observable instead of the fundamental constituents, quarks and gluons, which are 
confined. Ten years before the advent of QCD, Skyrme proposed a model [p]] for hadronic 
physics which involves only meson fields (pions), and where baryons emerge as topological 
solitons. This model is now recognized as the simplest chiral realization of QCD at low energy 
and large iV c 0. In the following we will refer to this model as the minimal Skyrme model. 

The minimal Skyrme model has been studied extensively to describe static properties of 
baryons [3-7] as well as dynamical properties, in particular the simplest vibrational excitation 
of the Skyrmion : the breathing mode [8-15]. Concerning the static properties, the predictions 
of the model are generally within 30% of experimental values when the parameters are adjusted 
to fit the nucleon and the A masses. In this paper we investigate whether this discrepancy is 
due to the use of the minimal Skyrme model. One motivation for this study is that it is easy 
to show that the second and the fourth order terms in the pion field derivatives have the same 
contribution to the soliton mass so that further terms should be added. Moreover, in taking 
into account the Casimir energy of the Skyrmion, the authors of Refs. [|I1| [TTJ have found that 



the nucleon mass is lowered to a value between and 400 MeV which is much too small ! They 
also proved that the O(N^) and 0(N®) mass contributions are of same order. Thus the large 
N c expansion is a poor approximation within this minimal Skyrme model. 

In this article we mainly focus on the breathing mode. One problem regarding this mode 
in the minimal Skyrme model is that it does not show up in analysis based on phase shifts 
i~0l [T^| , whereas when other calculations found it || |, [D], [13], [14], |15| , its excitation energy 



stands between 200 and 300 MeV which is too small compared to its experimental value of 500 
MeV. 

One possible way to circumvent these difficulties is to improve this model by adding higher 
order terms in derivatives in the corresponding chiral Lagrangian as the sixth order term gen- 
erated by u- meson exchange [jl7j (first proposed in Ref. [|i~8|j ). 

This extended model has already been used to describe the Roper resonance. Kaulfuss and 
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Meissner |I9| used both scaling approximation and semiclassical quantization. They found a 
resonance at an excitation energy of ~ 480 MeV. However, these authors used values for the 
parameters of the model which are in conflict with those determined by chiral perturbation 



theory [2C, 21 1 and, moreover, changed the sign of the so-called Skyrme term. Schwesinger and 
Weigel |EJ have also used phase shifts analysis within this extended Skyrme model. However, 
they did not find the Roper resonance. 

In this work, we describe this low-lying monopole resonance within the same extended 



Skyrme model of Refs. JT7 , |22| using the linear response theory. This method is more transpar- 
ent and has already been shown to be powerful in different domains such as giant resonances 
in nuclear physics |23| or nucleon polarizabilities in hedgehog models in hadronic physic |24[ . 
In a previous article |T5j we demonstrated the practicability of this approach. However this 
calculation was limited to the minimal Skyrme model only. 

The present article is organized as follows. In Sec. I we introduce the extended Skyrme 
model and define our notation. In Sec. II we review the linear response approach of Ref. [|TI| 
and specialize to the extended model considered here. We finally discuss our results concerning 
some static properties of the soliton and the Roper resonance in Sec. III. 



1 The extended Skyrme model 



The simplest chiral effective Lagrangian proposed by Skyrme [H reads 



C 



SK 
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Tr {d^Ud^W 



32e 2 



Tr{ UdJJ)U + ,{d v U)U + ] } 



where U is an SU(2) matrix parametrized by the (Goldstone) pion fields 7r a , normalized to the 
pion decay constant F n : 

U = exp[ 2m.r/F n ] , (1.2) 

the r a 's being the usual Pauli matrices. The first term in the Lagrangian ( |1.1| ) corresponds 
to the well known nonlinear a model and the second one, which is of fourth-order in powers 
of the derivatives of the pion field, was introduced by Skyrme to stabilize the soliton. It is 
generally referred to as the Skyrme term. It can be derived from a local approximation of 
an effective model with p-mesons |[25|| . Similarly, a term of order six can be generated from 
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u;-meson exchange [18|, |26| . It reads : 



1 3 2 

A = -^B^ 1.3 



where the anomalous baryon current is given by 

B, = ^^Tr { (d v U)U + (d a U)U + (dpU)U+ }. (1.4) 

The two new constants appearing in Eq. (|1.3| ) are the u;-meson mass and the parameter (3^ 
which can be related to the uo — > 7T7 width |p7fl . 

In Ref. |17[] it was shown that the effect of the Casimir energy, including the term Cq in the 
Lagrangian (p..l|) , is to lower the nucleon mass by about 500 MeV, so that one obtains a more 
satisfactory value of ~1 GeV. Anticipating on the last section, let us observe that the axial 
vector coupling constant g& is found to be 1.24 instead of the Skyrme model prediction 0.34. 
So there are good reasons to think that the extension of the Skyrme model which consist in 
adding the sixth order term ( |1.3j ) to the Lagrangian ( |1 . 1|) is a good approximation to describe 
the baryonic sector. 

There are of course other sixth order terms in the pion field derivatives such as the contri- 
butions of the p and scalar mesons. Nevertheless, a Lagrangian containing all possible sixth 
order terms leads in general to an Euler-Lagrange differential equation of order higher than 
two. Consequently, one is not sure to find a soliton type solution. This question is still open up 
to now since the parameters which correspond to these terms have not yet been determined. 
Fortunately, the term Cq in Eq. ( |1.3|) has the noteworthy property (as the Skyrme term) of 
being quadratic in the time derivatives and to lead to an equation of motion of second order as 
we will see below (see Eq. ( |1.7| )). For this reason we consider only this term in this work. 

Our starting point is then the following Lagrangian : 

C = C S K + C 6 + ^mlTr(u + U + -2) . (1.5) 

The first and second term in Eq. (|1.5| ) have been discussed above. The last term which is 
proportional to the square of the pion mass m n implements a small explicit breaking of chiral 
symmetry. Assuming the hedgehog ansatz for the pion field v?(r,t) = F n F(r,t) r/2 , (see Eq. 
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( |1.2| )) the Lagrangian density ( |1.5| ) becomes 



+ 



sin 2 F 



e 2 r 2 
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01 



/sin F 



sin F 
2r 



3m 2 7r 4 V r 
sin 2 F 



d^Fd^F 



g2 ,y«2 



+ -F>*(cosF-l) 



The corresponding classical Euler-Lagrange equation reads 



(1.6) 



g*F + aF* = ig 2 F' 



9 + aF 
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where the time-dependent functions g, a and 6 are, respectively, 



g(r,t) 

a(r, t) 

e(r,t). 
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r + sin F 



sin(2F) 1 . „ , 9 9 9 „ 
— + -(ei^m 2 r 2 sin F . 



;i.7) 



;i.8) 



Primes and dots in Eq. (|1.7|) indicate radial coordinate differentiations and time differentiations 
respectively. In order to ensure that the baryon number is equal to one, the equation of motion 
(11.71) has to be solved with the following conditions : 



F(0,t) 
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F(oo,t) = 



;i.9) 



which are sufficient for a differential equation of second order. 



2 Linear response analysis 

In order to describe the low-lying monopole vibrations of the Skyrmion within the extended 
Skyrme model (|1.5| ) we explore the response of the Skyrmion to an external infinitesimal 
monopole field with a frequency Q. The response function is determined from the evolution of 
the isoscalar mean square radius of the Skyrmion with respect to the frequency Q. 

An external time-dependent monopole field corresponds to the addition of the following 
term to the Lagrangian density Q1.6Q 



Ant = eF^r 2 B°(r,t) e sin (fit) exp(^t) 



(2.1) 
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where B°(r,t) is the time-component of the baryon current ( O ) 

and 77 a vanishingly small positive number. Adding this interaction term (|2.1| ) to the Lagrangian 
density ( |1.6|) , the new corresponding Euler-Lagrange equation reads 

g 2 F + aF 2 = (g 2 F')' - (o + aF' 2 ) + e^^- sin 2 Fsin(fit) exp(^) (2.2) 

where g, a and # have been defined in the previous section. This last equation is to be solved 
with the boundary conditions F(t = — oo, r) = F s (r) and F(t = — oo, r) = where F s (r) is the 
static solution of Eq. (|1.7| ). 

Because the term ( |2.1| ) is weak (in the domain t G [—oo, 0] ), it introduces small changes of 
the classical solution. We thus can treat this solution in a linear approximation. To first order 
in e, F(r, t) has the form 

F(r, t) = F s (r) + 5F(r, t) + 5F*(r, t) , 
where SF(r, t) is linear in the field strength e. Moreover, its time-dependence reads 

SF(r, t) = -i^R(r) exp(z(fi - irj)t) 



Thus, equation (|2.2|) becomes 



W-i v y-A s )(g s R)= ( e ^) 3rsin2 ^ 



7T 



(2.3) 



where the function g s and the operator A s are, respectively, 
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- - sin 2 F s 



m 2 r 2 cos F s 



(2.5) 
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The isoscalar mean square radius || is given by 

2 r°° 

(r 2 ) = / r 2 sin 2 (F)F' dr. 

7T JO 

Up to first order in e it reads 

(r 2 ) = (r 2 ) s - z|( f(n) exp(*(fi - irj)t) - /* (fi) exp(-i(Q + i V )t) ) 

where (r 2 ) s is the static mean square radius and / the response function 

4 r°° sin 2 F 1 

/(«) = -/ r ^)dr. (2.6) 

By using equation (|2.3|) , we can extract the following spectral representation of the response 
function 



/(Q) = --E 7?r^r J 2 > ( 2 - 7 ) 

where the state 4> is defined by 

/ u\ 2eF 7r rsin 2 F s 

( r 0) = , (2.8) 

and the n are the eigenstates of the operator A s (see Eq. (|2.5|)), with the eigenvalues u 2 , 
normalized according to 

/•oo 

4> n (r)(j)m(r)eF n dr = 5 nm . 



In Eq. ( |2.7| ) the limit 77 — ^ + is, as usual P3[ , implicit and corresponds to the boundary 
condition specified above. The quantity of interest here is the imaginary part of the response 
function, which is directly related to the distribution of collective strength (see e.g. [23| ). The 
energy at which the imaginary part of the response function (|2.7Q exhibits an unbound peak is 
identified with the excitation energy of the Roper resonance. 



3 Results and summary 

In this section we present our results for some static properties and for the energy of the Roper 
resonance within the extended Skyrme model ( |1.5| ) and compare them to those of the minimal 
one. First of all, because it is a meson theory, the parameters of the model have to be fixed by 
fitting the low energy meson observables. Concerning the pion decay constant, the pion and the 
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cu-mesons masses, the experiment yields 186 MeV, 139.5 MeV and 782 MeV, respectively. For 



the dimensionless parameter e, Riggenbach et al in analysing the Ku decays, reduce by a 
factor two the error bars on the chiral low-energy constants ||20|| . As a consequence, one finds 



e = 7.1 ±1.2 [H] which is surprisingly very close to the value e = 2ir proposed by Skyrme a long 
time ago [1|]. The last parameter (3^ is obtained by fitting the u — > ttj width, yielding /3 W = 9.3 
2"7 |. In order to compare the minimal and the extended Skyrme models we will consider two 
sets of parameters : 

[i] F n = 186 MeV , e = 7.1 and f3 u = (minimal Skyrme model ). 

[ii] F n = 186 MeV , e = 7.1 and (3^ = 9.3 (extended Skyrme model ). 

In Tab. I we report the soliton mass, the isoscalar root mean square radius, the axial vector 
coupling constant and the A-nucleon mass splitting for the different sets of parameters. In the 



third row we report the experimental values [28 



Concerning the soliton mass we find 972 MeV and 1563 MeV in the case [i] and [ii], re- 
spectively. We have to subtract ~ 1 GeV (Casimir energy) in the case of the minimal model 
(case [i]) which leads to a value of ~ MeV [T7J ! In the extended Skyrme model (case [ii]) 
we must subtract ~ 500 MeV 17] which leads to a reasonable value of ~ 1 GeV. Regarding 
the A-nucleon mass splitting, which is not affected by the Casimir effect, we find 1192 MeV in 
the case [i] which is not a realistic value compared to the experimental one 290 MeV. However, 
in the case [ii] we obtain 227 MeV which is more acceptable. Furthermore, the axial vector 
coupling constant gA is found to be 0.34, in the case [i], instead of the experimental value 1.23 
which is very close to the value predicted by the improved Skyrme model (case [ii]). Of course, 
one has to account for the Casimir effect in the calculation of gA- However, the corrections seem 



to be small within this extended Skyrme model JT7|]. For reference we plot in Fig. 1 the chiral 
function F s solution of the static Euler-Lagrange equation ( |1.7| ) for the minimal and extended 
model. 

For the breathing mode of the Skyrmion, we display in Fig. 2 the imaginary part of the 
response function ( [2.7|) for the two sets of parameters. In both cases we find an unbound peak. 
In case [i] we find this peak to be broad and located at 500 MeV. Nevertheless, this result 
cannot be reliable since the other predicted properties are unrealistic (see Tab.l). Contrary, 
in the case of the extended Skyrme model, the response function exhibits a pronounced sharp 
peak located at ~ 400 MeV which we identify to the Roper resonance. This value corresponds 
to the Roper- nucleon mass splitting and, consequently, is not sensitive to the Casimir effect. 
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As mentioned in the introduction, the authors of Ref. []22] found no trace of the Roper 
resonance within the same extended Skyrme model. So it is natural to check if this result is 
due to the parameters they chose. When we take their parameters ( = 142.4 MeV, e = 9.92 
and /3 W = 13.6 ) to investigate the Roper resonance with our approach, we find a sharp peak 
located at 320 MeV. This leads us to think that the difficulty in finding the breathing mode 
resonance in Ref. [2"2~ is due to the implementation of the phase shifts method. 

To summarize the main results of our calculations, we plot in Fig. 3 the energy spectrum of 
the A and Roper (N*) resonances for the different values of the parameters and compare them 
to the experimental one. We see obviously that the case [ii], which corresponds to the extended 
model, is the closest to the experimental situation. In order to obtain a better agreement one 
has to take into account higher order terms in addition to (see Eq. ( |L3p ) in the chiral 
Lagrangian. 

Finally, the message that we want to transmit through this work is that one should not 
restrict oneself to the standard Skyrme model |I] (see Eq. ( |1.1| ) ) for the description of low 
energy hadron physics, but consider extensions of this model including higher order terms 
in powers of the derivatives of the pion field. In this sense, the model considered here can 
be considered as a minimal extension of the Skyrme Lagrangian. This claim confirms the 



conclusions of Ref. |17|, [29| . A more realistic improvement consists of considering effective 
Lagrangians which incorporate low mass mesons with finite mass 
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Table captions 



TABLE I. Are reported, for the two different combinations of the parameters, the soliton 
mass Ms, the isoscalar root-mean-square radius r , the axial vector coupling constant g^, the 
A-nucleon mass splitting and the Roper-nucleon mass splitting hQR oper . The last line is derived 
from the data |28[| . F n , m n and are taken to their experimental values (see Sec. III). 



e 




M s 


r 


9a 


M A - M^v 








MeV 


fm 




MeV 


MeV 


7.lW 





972 
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0.72 


1.23 
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a) from Ref. pi and (b) from Ref. [p7f 



Figure captions 



FIG. 1. Static hedgehog solution F s (r) with experimental values of F n , m n and m u (see 
Sec. III). 

FIG. 2. Imaginary part of the response function / (fm 2 ) versus the energy Q in MeV with 
experimental values of F n , m w and m u (see Sec. III). 

FIG. 3. Energies splitting of the A and Roper iV* resonances with respect to the nucleon 



iV according to the cases [i] and [ii] of the parameters (see Sec. Ill) and to the experiment [28 
All the energies are in MeV. 



